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ABSTRACT 


A simplified hydrocarbon-air combustion gas model, including the 
effects of dissociation, for convenient use in engine cycle computer 
programs is presented. The gene~alized model reduces to the hydrogen- 
air system as well as dissociating air. The exclusion of chemical 
species containing atomic nitrogen allows a considerable simplification 
of the composition equations. The thermodynamic properties of stoichio- 
metric combustion of the kerosene-air and hydrogen-air systems are 
computed with the simplified model and compared with more comprehensive 
gas models. In addition, the effect of the neglected chemicals species 
on the performance of a idealized subsonic combustion ramjet is presented. 
The simplified gas model has been used to define the limiting conditions 
for solid carbon and ammonia formation for fuel rich gas mixtures. A 
computer program listing of the simplified gas model calculation 
procedure is presented. 
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acceleration due to gravitational field of the earth, 
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ratio of hydrogen atom to carbon atom in fuel molecule 
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Subscripts 


» free -stream conditions 

J Jth species 

J Jet 

lim limiting condition for solid carbon or ammonia formation 

s static conditions 

t stagnation conditions 

Superscripts 

0 at reference pressure, 1 atmosphere 


CHAPTER I 


INTRODUCTION 

Recent and continuing advances in fluid dynamics and metallurgy 
have opened new domains of power and efficiency to the air-breathing 
combustion engine. The combined effect of technological advances and 
the demand for greater operating flexibility in the aircraft power 
plants have led to certain design trends which are significant from a 
thermodynamic standpoint, namely: the maximum allowable internal 

temperatures, pressures, and gas stream velocities are increasing, and 
minimum internal operating pressures sure decreasing as well. 

Developments in the area of compressor and turbine blade cooling 
have been instrumental in extending the potential range of efficient 
operation of the turbojet engine to high supersonic Mach numbers and 
high altitudes. As discussed in Zucrow (Ref. l), for turbojets 
employing uncooled turbine blades, the temperature of the gases leaving 
the combustion chambers must be limited to approximately 2000° R due to 
stress, creep, and stress duration considerations. Because of that 
temperature limitation the overall fuel-air ratio for a turbojet engine 
burning a hydrocarbon fuel, such as JP-k (CH 2 ), is quite low, of the 
order of 0.015 (R eq *» 0.25). Current turbojet designs, such as those 
proposed for the supersonic transport, employ turbine blade cooling 
which permits turbine inlet temperatures of approximately 2800 r R 
(Rothrock, Ref. 2) and results in fuel-air ratios of the order of 0.03 
( R eq ** Advanced turbojet engines which employ hydrogen for both 
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fuel and turbine cooling are being studied with stoichiometric combustion 
and resulting turbine inlet temperatures of approximately 4000° R. 

The ramjet engine is better suited for propulsion problems 
requiring high flight speeds. No such temperature limitation is 
applicable to the ramjet engine because it has no moving parts that are 
subjected simultaneously to dynamic loads and high temperatures. Con- 
sequently, stoichiometric combustion of high energy fuels, such as 
methane (CH^) and hydrogen, may be utilized for the ramjet engine. In 
addition, the high stagnation temperatures attained by the air at high 
flight Mach numbers result in combustion temperatures approaching 
5000° R for stoichiometric methane combustion and 6500° R for stoichio- 
metric hydrogen combustion. 

As discussed by Ferri (Ref. 3), the concept of the supersonic 
combustion ramjet appears attractive for very high Mach numbers and 
altitudes approaching orbital conditions. Since the fuel is mixed and 
ignited in a supersonic airstream, internal temperatures would not 
exceed 5000° R (stoichiometric hydrogen-air combustion), even at flight 
Mach numbers of the order of 12. Because of the application of the 
supersonic combustion ramjet to high altitude flight and the fact that 
a large portion of the energy in the air6tream remains in the form of 
kinetic energy, internal operating pressures can be substantially below 
sea level atmospheric pressure. 

An accurate knowledge of the thermodynomic properties and 
behavior of these engine’s working fluid under all operating conditions 
is a basic requirement for effective design. However, at high temperatures 


3 


and low pressures, considerable difficulties arise in calculating the 
thermodynamic properties of combustion gas mixtures - difficulties which 
stem from the phenomena of internal energy excitation coupled with 
chemical dissociation. 

In the past, methods such as presented by Hall and Weber in 
Reference 4 have been used to calculate the thermodynamic properties of 
combustion gases for air-breathing engine cycle computations. The method 
of Reference 4 is based on separate calculations of thermodynamic 
properties of air and the products of combustion of a stoichiometric 
fuel-air mixture . It is assumed that the properties of the combustion 
products for any fuel-air ratio less than stoichiometric may be obtained 
by linear interpolation between the two extreme cases. Results of this 
method are exact for cases with no dissociation. However, when disso- 
ciation is taken into account, errors arise because the method does not 
account for the change in composition of the gases when they axe mixed. 
This error is small for turbojet or turbofan calculations since tempera- 
tures are relatively low a- * dissociation is not extensive. However, 
the high operating temperatures of hypersonic engine cycles may result 
in cases of extreme dissociation and cause large errors in calculating 
engine performance. 

The method of Reference 4 does not allow engine calculations for 
equivalence ratios greater than unity. However, at high Mach number 
flight, the fuel required to cool aerodynamic surfaces and engine 
component may force engine operation into this region. 
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In short, the method of Reference ^ does -'vt provide a sufficiently 
general basis for engine cycle calculations under all conditions of 
current and future interest. 

Chemical equilibrium theory which governs dissociation phenomena 
has long been formulated, for example by Clarke and McChesney 
(Ref. 5 )* However, its application to complex combustion gas mixtures 
has proven so tedious and complicated as to preclude direct use in 
engine cycle computer programs. The General Electric combustion gas 
models of References 6 and J, which eerve as a basis of comparison for 
the simplified gas model presented herein, are elaborate treatments of 
the hydrogen-air and kerosene -air systems . 

The generalized hydrogen -air model of Reference 6 assumes 
12 chemical species,, requiring eight independent chemical reactions 
leading to eight equilibrium expressions. Since the initial proportions 
of hydrogen and air define four mass balance equations, the mathematic. 1 
solution: 12 equations - 12 unknowns, is demonstrated. The reduction 

of the system of equations to one equation - one unknown, as recommended 
by Erickson, Kemper, and Allison in Reference 8 - is extremely tedious. 
With the exclusion of the technique of Reference 8, hope for the 
solution of this system by a single level iteration is abandoned. The 
adoption of a two-level iteration or one of the methods summarized by 
Zeleznik and Gordon in Reference 9 to solve for gas composition has a 
great effect on the utility of the computation procedure as a subroutine 
for an engine cycle computer program. Consequently, the most efficient 
means of representing the thermodynamic properties of these gas models 
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in an engine cycle program is by an elaborate system of fitted curves, 
sacrificing both computer storage and accuracy. 

The purpose of the simplified gas model presented herein is to 
incorporate the substantial effects of dissociation consistent with 
convenient use in computerized engine cycle calculations. The proposed 
gas model is simplified in the sense that the formation of species 
containing atomic nitrogen is neglected. This assumption enables the 
solution for chemical composition to be obtained with a single level 
iteration. Although the original intent of this study was directed 
toward the hydrogen-air system, it was found that the model could be 
generalized to any hydrocarbon -air system, as well as dissociating air, 
with little additional complexity. The computer program included can 
be readily incorporated as a subroutine in an engine *ycle program or 
used alone to generate Mollier diagrams . Inputs to the program are the 
ratio of carbon atoms to hydrogen atoms in the fuel molecule, equivalence 
ratio, temperature, and pressure. Outputs from the program are the mole 
fractions of the chemical species assumed, enthalpy, entropy, and 
molecular weight. 


CHAPTER II 


ANALYSIS 

Description of Gas Model 

The range of temperatures and pressures of interest for air- 
breathing engine application is assumed to be below temperatures of 
7000° R and between pressures of 0.001 and 100 atmospheres. In this 
range of conditions. References 3 and ^ indicate that the formation of 
nitrogen species, such as N, NH, NH^, and NO occurs in negligible 
amounts and, therefore, has a very small effect on the thermodynamic 
properties of combustion gas mixtures. The assumption to neglect these 
species has a great effect on the complexity of the calculation 
procedure since molecular nitrogen can be considered inert. The effect 
of this assumption on the thermodynamic properties of gas mixtures and 
on the performance of an idealized ramjet engine will be discussed in 
the section entitled Results and Discussion. 

Figure 1 is a diagram illustrating the chemical species and 

required chemical reactions considered in this gas model. The figure 

is divided into regions of temperature and equivalence ratio. Above 

T . „ the gas is considered to be dissociating. The chemical 

cut-off 

species assumed for pure air (R e q = 0) are identical to those of 
Reference 6 and include atomic nitrogen species. For the combustion 
gas model (R > 0) the dissociated nitrogen species (N, NH, NHj, NO) 
are not included. It is apparent that by excluding the species and 
reaction containing carbon, the model reduces to the hydrogen-air system. 


■sm 
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Figure 1.- Description of gas model. 
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Below T? cu £_ 0 ff, in general, all chemical reactions are complete 
and combustion products completely formed. Therefore, the initial 
proportions of elements are sufficient to define the gas composition. 
However, for R e ^ > 1, the gas model of Reference 7 indicates that CO 
may form even below T cut _ 0 f f . As a result, the initial proportions of 
elements are not sufficient to define the gas composition and one chemical 
reaction is required. 

In the region of fuel rich operation, solid carbon and ammonia 
may form. The limiting pressure for the formation of solid carbon 
depends upon the relative proportions of carbon dioxide (CO 2 ) and carbon 
monoxide (CO), equivalence ratio, and temperature. However, there is a 
limiting equivalence ratio for a given hydrocarbon-air system above 
which carbon will form under any condition of pressure and temperature. 

The formation of ammonia was considered only for the hydrogen-air 
case since the formation of solid carbon preceded the formation of 
ammonia for the combustion products of all the hydrocarbon-air systems 
studied. The limiting pressure for the formation of a given amount of 
ammonia depends upon the relative proportion of nitrogen (N 2 ) and 
hydrogen (E^), equivalence ratio, and temperature. 

It should be emphasized that the combustion gas model, including 
the criteria for solid ca-bon and ammonia formation, is based entirely 
upon equilibrium theory of gases. 

Derivation of Gas Model 

The general equations presented in this section for computing 
chemical composition and combustion gas properties make use of the 
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following assumptions: (l) the formulated gas mixtures simulate the 
products of complete, constant pressure, adiabatic combustion; (2) the 
gas mixtures are in thermal and chemical equilibrium; (3) the ideal gas 
equation of state is valid; (4) only gas phases are considered; and 
(5) the effects of ionization sure negligible. 

Overall stoichiometry and gas model constants .- In order that the 
simplified gas model may apply to the combustion of any hydrocarbon fuel 
a generalized statement concerning the atomic composition of the fuel is 
required. Consider the following stoichiometric reactions: 

CH + ^ 0 2 = C0 2 + | H 2 0 

CH 2 + f °2 = c0 2 + | h 2° 

CH, + l 0 2 = C0 2 + | 11,0 

CH4 + | o 2 — co 2 + | H 2 0 

In general terms: 

CBn + 0 2 = C0 2 + | H 2 0 

where 

n = no ‘ atoms H in fuel molecule 
no. atoms C 

Thus, in the general case; 

no. moles CHq 4R e q 

no. moles 0 2 n + 4 
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where R e q Is the equivalence ratio. The initial composition of air 
incorporated in this gas model is the same as that of Reference 6; that 
is, combustion is assumed to have occurred with dry air of the following 
composition by volume: 

0 2 N 2 A 

20 . 9495 per cent 78.0881 per cent 0.9624 per cent 

The following constants and parameters define, in general terms, the 
initial atomic proportions of the elements in the gas mixture. 


c = no. atoms A = 0 . 0 06l6227 
x no. atoms N 


. no. atoms K , 5.727445 
2 no. atoms 0 


p _ no. atoms C _ ^ e q 
^ no. atoms 0 n + 4 


c 4 “ 


no. atoms H _ ^ n ^ e< l 
no • at- ops 0 n + 4 


Defining R^ = i, then 


C - 

^ 1 + 4RQJJ 


= 


2Req 

1 + ^h 
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For the special case of hydrogen-air combustion: 


2H2 + O2 " 1 2H2O 

no. moles Hg _ 2R 
no. moles O 2 e< * 

Rch = 0 Cj = 0 C = 2R e q 

The stoichiometric fuel-air ratio is given by: 


f s = 0.028931 


1.008 + 12 . 031 ^ 


1 + bR 


ch 


It is of value to represent the equilibrium constants for the 
assumed reactions in a convenient manner. The treatment of chemical 
equilibria in Reference 5 leads to the equilibrium expression for a 
general reaction: 


aA + bB cC + dD 


The equilibrium expression is: 


Xc*D p-(c+d-a-b) . K 

a a a b 


where 


K = exp 



a l 

RT 




- b 
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Figure 2 shows the variation of log^g K with ^ for the chemical 

reactions considered in this gas model. The values of Gibbs free energy 
G° 

-ji for the species were taken from the data tabulations in References 6 
RT 

and 7* Figure 2 indicates that the equilibrium constants can be 
represented as linear functions. 

Hydrocarbon-air gas model with dissociation .- The following gas 
model applies to the region: 


T > T 


cut-off* 



> 0 


Chemical species assumed: 


HgO, CC>2 y A, Oj Hj CO j OH 


Chemical reactions: 


2H + 0 = EgO, 2H = Hg, 20 = 0 2 , 0 + CO = C0 2 , H + 0 = OH 


Mass conservation equations: 

10 

Z - 1 

J=1 


( 1 ) 


X A = C 1 (2X N2 ) 


( 2 ) 


^ - C 2 (2Xco 2 + *00 + V + ^ + *0 + X 0 H ) 


( 5 ) 


lOOr 


80 h 


60 r~ 


o 40- 

o< 

o 

20 - 


Oh 


-20 L- 


Figure 2. 


•13 



0 .0002 .0004 .0006 .0008 .001 

l/T, l/°R 


i i i i 

7000 3000 2000 1000 

Temperoture,°R 

3000 2000 1000 600 

Temperoture,°K 


- Variation of equilibrium constants with temperature for gas 
model reactions. 
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Xcc>2 + ^0 = C 3( 2X C0 2 + X C0 + x H 2 0 + ^Og + x 0 + X 0H^ 


X H + ^2 + 2X H20 + X 0H = C 4^ 2X C02 + X C0 + X H20 + X 0 + X 0H + 2X 0 2 ^ 


(5) 


Equilibrium expressions: 


x h 2 o - ^Wh (6) 

x Hg = (T) 

Xq 2 = PKjXg (8) 


x C0g = PX 4 X O X C0 


(9) 


Xq H = PK^XqXjj 


( 10 ) 


Solving for two equations - two unknowns in X Q , X^: 


a * {a 2 2 X 0 + ^^0 + A20} X i + { a 12 X § + A 11 X 0 + A 1q) X H 
+ + *02$ * Aoi 3 ^ + fk oo} a 0 
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p = (b 22 X§ + B 21 X 0 + B 2 q)x| + (b 12 X§ + B 11 X 0 + 3 10 }x h 
- {BojX^ + 1^2*0 + ^)lXo> = o 

For R g(1 > 1 the iteration function, a, is insensitive to 
changes in Xq. As a result, Xg is used as the iteration variable in 
this region of computation. Unfortunately, if Xg is assumed, the 
expression p is cubic in Xq and an additional iteration is required. 
For fuel 'rich hydrogen -air mixtures =0, p is quadratic in Xq, 

and the single level iteration is maintained. However, for fuel rich 
mixtures containing carbon, a double level iteration is necessary for 
solution, where 

h.2 “ p5k iM1 - 2C 5 + C) 

Agi = P^KjU + C) + 2X^(1 - 2C 5 )} 

A 2 q = 2PK 2 (l - Cj ) 

A 12 = P^K^l - 2C 5 + C) 

A 11 = KV 1 + c ) + 2*4 (1 - 2C 3 )} 

A 10 = " C 3^ 

Aoj = 2CP 2 K 3 K4 

Aq 2 = P{2K 3 (C + Cj) + K^(l - 2Cj + C)) 
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Aoi = 1 + C - 2PK 4 (l - 2C ? ) 

A 00 = " 2 ( 1 " 

and 

B 22 = P3K 1 K 4 {2(1 - C 5 ) - C 4 ) 

B 2i = P 2 {k x | 2(1 - C 3 ) - Cj + 2X^(1 - 2C 5 )} 

B 20 = 2PK2(1 - C 3 ) 

B 12 = PV^l - 2C 5 - 0 k ) 

B ll = p ( K 5 (1 " C 3 ~ + V 1 " 2C 3 ] ) 

B 10 = 1 - c 3 

^3 = 2C 4 P% 3 K 4 

Bq 2 = C 4 P{2K 5 + K4) 

BOl = C 4 

C = 1 + C 2 (l + 2C^) 

Setting C 3 and equal to zero reduces the system to identically 

the hydrogen-air case. The expressions for a and 0 are satisfied 
simultaneously by use of the Newton Iteration Scheme. Newton’s iteration 
scheme for determining X Q and Xg can be written as 
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X 0 (L + 1) - X 0 (L) - 


g(Xo(L),X H (L)) 

~(a(X 0 (L),X H (L))) 


Where Xq(L) is the Lth approximation of Xq and the derivative of 
a(X 0 , Xg) with respect to Xq is 


da da bX Q 
dX Q " dX H dp 
dX H 


Air-gas model with dissociation.- The following gas model applies 


to the region: 


^ ^cut-off* ^eq ® 


Chemical species assumed: 


0 2 , N 2 , A, 0, NO, N 


Reactions: 


20 = 0 2 , 2N = N 2 , N + 0 = NO 


Mass conservation equations 


X J = 1 


( 1 ) 
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X A * c l( 2X N2 + X N0 + X N^ ( 2 ) 

2Xn 2 + X N0 + X N = C 2( 2X 0 2 + X 0 + X NC>) (5) 

Equilibrium expressions 

^ w 

% 2 - ( 5 ) 

X N0 = 2 X 7 X 0 % ( 6 ) 

Solving for two equations two unknowns in Xq, Xjj: 


a = + { A iiXo + a io) x n + fyo 2 x o + ^ 1 % + A oo} " 0 

0 = B 2Q x| + {b-qXq + B 10 }X N - {Bo 2 Xq + ^)l x 0> = 0 

where 


A 20 = (1 + 2C!)PK6 

b 20 = 2PK 6 

A U = (i + C^PKr 

®11 “ (1 “ C2)PKy 

A 10 “ 1 + C 1 

H W 

o 

11 

H 

A02 * PIC 3 

%> “ 2C2PK, 

AO!* 1 

Boi = C 2 

AOO - -1 
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The simultaneous solution of a and 0 is subject to the same iteration 
procedure as the gas model of the previous section. 

Fuel rich gas model with solid carbon formation .- The following 
gas model applies to the region: 

T - T cut-off> R eq > 1 

Chemical species assumed: 

1^0, CO gj Hgj N 2 j A f CO 

Chemical reactions: 


Eq + COg — CO + HgO 
Mass conservation equations: 


6 



J=1 



( 2 ) 


“ C 2^ 2X C0 2 + X C0 + X HoO) (5) 

00 


X C0 2 + X C0 = C 3^ 2X C0 2 + x co + x h 2 o) 
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2XH 2 + 2XH 2 0 - C4(2Xco 2 + X co + X^q) (5) 

Equilibrium expressions: 

V^O - %Vc°2 (6) 

Solving for Xqq yields: 

AAX§ 0 + BHX C0 + CC = 0 

where 

AA = (1 - Kq)(C + + 2Cj - l) 2 

BB = (C + C k + 2Cj - 1) (2(1 - 2C 5 ) + Kq(C4 + 6C ? - 2)J 
CC = 2CjKq(2 - Cj,. - kCj) 

Consider the following reaction: 

C0 2 + C(S) = 2C0 

Equilibrium expression: 


P lim = 


K^COg 

*80 


If P > P lim , solid carbon will form. There is, however, an equivalence 
ratio for the combustion products of a given hydrocarbon (R c j 1 ) , where 
solid carbon will form even in a vacuum. 
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Consider the equilibrium expression: 


P lim X CO " K 9 X C0 2 


p Hm X C0~ 


Setting = 0 in the system last derived, gives: 


3 _ 1 + ^ R ch 

eq, lim 2 R ch 


Nondissociating gas model.- The following gas model applies to 


the region: 


T 1 T cut-off 


Chemical species assumed: 


HgO, CO 2 , 0 2 , N 2 , ^ 


Mass conservation equations: 


x j = 1 


X A - ^(2% 2 ) 


2%, * C2<2XC0 2 + X H20 ♦ 2Xq,) 
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= C 3^ 2X C02 + X H 2 0 + 2x 0 2 3 


2X H 2 + 2X H 2 0 ~ c 4( 2X C0 2 + X H 2 0 + 2X 0 ? ) 


For 1.0 > R eq > 0 


x H 2 » Xj - Xj(R ch ,R eq ) 


For Req «= 1 


X H 2 " °> ^ = °> X J = x j(R ch ) 


For the hydrogen-air case (R ch = 0) R eq > 1 


X 02 “ Xj = Xj(R e q) 


consider the following reaction: 


N 2 + 3H 2 = r W-. 


Equilibrium expression 


P H“ ‘ X S 2 X »2 


X BH 5 >0-01' P > P li„ 
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Thermodynamic properties of gas mixtures .- The previous deriva- 
tions have outlined the solution for gas mixture composition. With the 
chemical composition defined and a tabulation of the thermodynamic 
properties of the pure constituents provided, the thermodynamic properties 
of a gas mixture can be calculated by the following (Ref. 6): 



Mean Molecular Weight: M = ZjXjMj 



constituents at reference pressure (l atmosphere). The thermodynamic 
expressions above inherently assume that the pure constituents obey the 
perfect gas law. 

The standard reference state of the elements A, N, 0, H, C are 
taken from Reference 6, that is: 

A as A 
N as N 2 
0 as 
H as Hg 
C as C (solid) 


24 


By definition, the energy content of these elements in their standard 
reference states (T = 0° R) is zero. Equations for the computation of 
fuel enthalpies consistent with this thermodynamic basis are presented 
in Reference 7 . 

A listing of the computer program to calculate the composition 
and thermodynamic properties of this gas model is presented and discussed 
in the Appendix. Computational time for this program has been estimated 
at 6000 cases /min on the IBM 7094 data processing system at the Langley 
Research Center. One case is defined as the computation of one point 
on a Mollier Diagram. 



CHAPTER III 


RESULTS AND DISCUSSION 

Comparison of gas models .- The computer program for the simplified 
combustion gas model has been used to calculate the thermodynamic 
properties of the stoichiometric kerosene-air and hydrogen-air systems. 
The purpose of this calculation is to compare the results of the simpli- 
fied gas model with the more extensive trectments of References 6 and 7 . 
Since the thermodynamic properties of the pure constituents used in the 
simplified model were taken from these references, the differences 
between the results of the simplified model and the reference models 
are due to the formation of the neglected chemical species N, NH, NH^, 

NO. 

Figure 3 is a Mollier diagram for the stoichiometric products of 
combustion of kerosene with air (Rcjj = 0.5)* The solid curves are 
values plotted from the tabulated data of Reference 7. The dashed 
curves are values obtained from the computer program for the simplified 
model. Unfortunately, the data presented in Reference 7 are limited to 
temperatures below 5000° R. Since the formation of species includi.g 
atomic nitrogen (N, NH, NHj) is more important at temperatures ab ve 
5000° R, the good agreement between gas models is not surprising. 

Figure 4 is a Mollier diagram for the stoichiometric products of 
combustion of hydrogen with air (R^ ® 0). The solid curves are values 
plotted from the tabulated data of Reference 6 . The dashed curves are 
values obtained from the computer program for the simplified model. 


25 



Entholpy, H, Btu/lbm 


26 



1 i 1 1 1 1 1 i . 

6 7 8 9 10 II 12 13X10* 

Eotropy.S, J/kg-»K 


Figure 3«- Mollier diagram; kerosene -air combustion products 

(R ch * 0.5), R eq = 1.0. 




Mollier diagram; hydrogen-air combustion product 
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The thermodynamic properties tabulated in Reference 6 cover temperatures 
as high as 10,000° R. The substantial disagreement above 7000° R at low 
pressures is due to the formation of atomic nitrogen species. An 
interesting result is that the agreement between lines of constant 
pressure is much better than the agreement between lines of constant 
temperature. This result is fortunate, since, for thermodynamic processes 
such as isentropic nozzle expansion, temperature is not used directly to 
calculate performance. 

The simultaneous conditions of high temperatures and low pressures, 
the area of substantial disagreement between the simplified gas model and 
the reference gas models, are beyond the realm of operation of typical 
hypersonic engine cycles. For example, for subsonic combustion ramjets, 
comristion pressures and temperatures are high since the airstream is 
brought nearly to stagnation conditions before combustion. For super- 
sonic combustion ramjets, combustion pressures and temperatures are low, 
since a large portion of the total air enthalpy remains in the form of 
kinetic energy. Hence, the conditions of high temperatures and low 
pressures do not occur simultaneously for engine cycles currently 
considered feasible for hypersonic flight. 

Ramjet performance comparison . - In applying the simplified gas 
model to hypersonic engine calculations, it is of interest to determine 
the effect which the use of the simplified gas model has on the cal- 
culated performance of a typical subsonic combustion ramjet. 

The subsonic c embus t ion ramjet is perhaps the simplest air- 
breathing engine cycle. Referring to the sketch below, the free -stream 
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air is decelerated through an air intake or inlet to low subsonic speeds 
in the burner. Here fuel is added at low velocity, combustion takes 
place, and then the flow is accelerated through a nozzle. 



Since the purpose of this performance calculation is to show the 
effect of small differences in thermodynamic properties on ramjet 
performance, the following simplifying assumptions were made: 

1. Free -stream Mach number chosen along a constant dynamic 
pressure path (q = 1500 psf). 

2. Airstream decelerated to stagnation conditions with total 
pressure recovery degraded to 10 atmospheres. (A good assumption at 
high Mach numbers, due to internal duct pressure limitations.) 

5. Completely mixed stoichiometric hydrogen-air combustion. 

4. No pressure losses during combustion (Pt nozz i e = 10 atmospheres). 

5. Enthalpy of injected hydrogen equals zero. Molecular hydrogen 

is a reference element, thus = 0 at T = 0° R. 

6. Combustion gas is isentropically expanded to free stream, 

' -utic conditions using so-celled "shifting equilibrium." 


30 


The 1962 standard atmosphere (Ref. 10) yields free-strearo static 
conditions, and together with assumption (l) gives the variation of 
altitude with Mach number shown in Figure 5* For steady adiabatic flow, 
the total energy in the airstream is given by: 

V§ 

H t,air ” H co,air + 2Jg 
For an adiabatic combustion process: 

^ _ ^air + ^^hydrogen 

^products 1 + f 


With two properties (H*. , , P + , ) of the combustion gas 

^products ^nozzle 

defined, the Mollier diagram can be entered. Figure 6 is a schematic of 
a Mollier diagram showing lines of constant pressure for the simplified 
gas model and that of Reference 6. Since the simplified gas model 
neglects the dissociated nitrogen species, the entropy level is slightly 
less than that obtained from Reference 6. In expanding the combustion 
gas to free-stream static pressure, the gas enters the region of almost 
exact agreement between gas models. Consequently, the expansion using 
the simplified gas model gives a slightly larger value of - H g than 
Reference 6. However, the effect of the difference in AH is reduced 
by the fact that: 
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APPENDIX 


COMPUTER PROGRAM 

The calculation procedure for determining gas composition and 
thermodynamic properties for the simplified gas model has been programed 
for and computed on the IBM 7094 data processing system at the Langley 
Research Center. A printout of the program is presented in this 
appendix. The program is written in Fortran language (Ref. 10). The 
sumbols used for the program sire as follows: 


RCH 

R ch 

X(l) 

x h 2 o 

REQ 

^eq 

X(2) 

x co 2 

TEMP 

T 

X(3) 

X °2 

P 

P 

X(4) 

Xh 2 

M 

M 

x(5) 

x n 2 

H 

H 

x(6) 

X A 

S 

S 

x(7) 

x 0 

MT(J) 

m j 

X(8) 

X H 

T7T(l) 

T i 

X(9) 

X C0 

HT(I,J) 

$ 

X(10) 

X 0H 

ST(I,J) 


x(ll) 

X N0 



X(12) 

X N 

Input.- The input Is 

read into 

the IBM 7094 

data processing 


system by the Fortran statement: 


1*6 


47 

READ (5,100) RCH, REQ, TEMP, P 
100 FORMAT (4E 16.8) 

For example, an input card would be: 

Col.- 1 17 25 49 

+0.00000000E+00 +0. 10000000E+01 +0. 50000000E+04 +0. 10000000E+00 

Output.- The output for this program is the mole fractions of 
each constituent and the thermodynamic properties of the gas mixture. 

For example, output for the case above would be: 

RCH=O.OOOOOOOOE-58 REQ=0.10000000E+01 TEMP=0.50000000E 04R P=0.100000000E 00 AIM 

XH20 XC02 X02 XH2 

O.I8515569E 00 O.OOOOOO0JE-38 0.27084661E-01 G.71114589E-01 

XN2 XA XO 

0.57613321L 00 0.71005768E-02 0.25695777E-01 0.63610586E-01 

XCO XOH XNO XN 

O.OOOOOOOOE-38 0.46109576E-01 O.OOOOOOOOE-38 O.OOOOOOOOE-38 

N=0.21993865E 02 H=0.l43l6l46E 04BTU/LBM S=0.309l8686E OIBTU/LBM- 

Built-in data.- The thermodynamic properties of the pure constit- 

f Hj s<?\ 

uents l—, — I are taken from References 6 and 7 and input as two- 
dimensional arrays named: 

HT(I,J) and ST(l,J) 


Where "I" is the index on temperature and "J" is the index on constit- 
uents. The thermodynamic properties are input for the discrete 


temperatures given in TT(l). The enthalpy array HT(l,J) is interpolated 
linearly while the entropy array ST(l, J) is interpolated logarithmically. 

Complete program .- The following program has been used on the 
IBM 709^ data processing system at the Langley Research Center. 




c thermodynamic properties of hydrocarbon-air combustion prooucts 

C H20 C02 02 H2 N2 A 0 H CO OH NO N 

C X ( I > X(2) XC3) X<4) X ( 5 > X < 6 1 XC7) X<8) X<9) XI101 XU1) XU2> 

DIMENSION TTU6 )*MT( 12) «HT ( 16f I2)«ST< |6i 12) *XU2 ) •Y«99) tF«99) • 

1 Z < 99 ) • TY ( 99 ) 

REAL M •’MT * K 1 #K2 *K3* KA *K5*K6 *K7*K8 * K9 *K l 0 
C 

C THERMODYNAMIC PROPERTIES OF PURE CONSTITUENTS 
C 

DATAIMT! J) * J=l * 12 1/ 1 8 • 0 1 6 • 44 *0 1 * 32 • 00 *2 *0 l 6 • 20 *0 1 6 • 

139*944 • 16*00* 1 *008* 28* 01* 1 7*000 *30 • 008 * 14*008/ 

DATA ( TT < t ) * 1 ■ l « 16 I/O* ♦ 200* *400* *600* *800*# 1000*t 1300* « 1600* #2000* ♦ 
12500* *3000* 1 3500* *4000* * 5000* * 6000 • « 7000 •/ 

DATA < ( HT ( I • J ) * 1*1 «16)«J*1 *5 1/-51727* *-51 100* • -50400* * -49700 • ♦ 

1- 48900* #-47743* #-46297* *-44939* * -42944 * « -40290* • -3744V# * -34466# ♦ 

2- 31354* *-24089* #-18176* #-11296* *-051 i 7. # -84000* # -03oOO* ♦ -82001 • * 

3- 81018* #-80 /46* #-79007. #-77147* #-74o3b* #-71 1 1 O* ♦ -67b9o* i-6JVcJj* • 

4- 6031 2* #-52057* *-44500* #-37000* *0* #698* # 1400* *2104* #2826* *35 /4* * 

54 744 * *5963* * 7639* #9794* # 1 1996* • 1424 1 • • 1 6524* *21203* *26021 • *30953* * 
60* *751 * « 1 356* #205<f* *2755* *3458* *4518* *5588* • 7042**0918* t 10867* • 

71 2086* •* 14964* * 192 72* *23744* #28355* *0* #699* • 1398* #2099* *2001 • •' 
83513**4608* *5743* #7316* *9356* • 1 1456* * 13597* * 

915768. *20171 .*24629. #29124*/ . 

DATA ( C HT (1*J)*!*1«16)« J*6*9)/0* *500* • 1000* *1500**2000* *2500* ♦ 


1 3250 • , 4090 • , 500C • « 6250 • , 7500. « 67b0 • * 1 3000* , 1 2300 • » 1 5000 • * 1 7500* « 

253425.. 34000., ?44J6., 54900. *553*5. , 56165. *56601 * .57600. . 55594. « 

« 

359655. . 61 11 0. *62356 . .63612. .661 20. *66640. . 71 J b3. *46759. *47300. . 

447712. . 48360 ».*fc'fe 70. ,4 9336. .49977. .50763. .51 752. *53006* * 54261 . . 
5555C3. .56756.. 69255. .61758. .64257. .-24639. .-24000. .-23000. . 

6- 2254 0.. -21 83?. .-21 life. *-20009. .-1 6656. .-1 7261 . .-1 51 93. . -1 3073. . 

7- 1 C91 2* ♦ -£726, * -4295. .0. *3500./ 

CATAJ (HTU.J >, 1*1, 16). J«1 0.12 J/8424. 9. 9098. .9750. .10455. .1 1222.. 

1 1 1 997 . . 1 004 3 . . 1 4 1 30 . ♦ 1 5605 . . 1 7520 . . 1 9503 .. 21 558 . * 23664 . . 27996 . • 

222456.. 3702a .. 19453. .20006. .20567. .21427. .22220. .23091 . .24248. * 

325445. . 2699*.. .29185 . ,31 341 • ,33580. .35743. .40220. .44741 . .49302. « 

4 1 01 022.. 102000.. 102000.. 102156. .1033 12.. 104590.. 105104*. 105933.. 
51 069C8. . I 08173. . 1 09422. .11 0663. . 1 1 1 916. . 1 1 4425. , 1 16946. • 1 1 951 8./ 

CATA < (ST 1 1 , J>« 1*1 .16). J«1 .4 I/O. * 18.68.21 .55*23.167,24.338,25.273. 
1 26. 434. 27. 40C.2e. 505,29 .691 ,30.725.31 .6464 ,32.4745.33 .91 63 .35. 1 >» 

236. 2003.0. . 21 .700. 24 .550. 26.20! 15.27.61677.28.61 060.30*32729. 

331 .61 326 .33. 0695 5. 34 .59469.35.5764 1 ,36.96959, 37.97133.39 .63861 » 

441 .0. 42. 19*0. ,21 .1536.23.6633.25.0721 .26.1033,26.9292.27.9544, 
528.7966.29.7306.30.6915,31 .4947,32.1666.32.796.33.8396.34.7178. 

635. 4778.0. .12. 5375. 14.7172. 16.1042.17.1073.17.8629. 18.8146.19.35. 
72C.3628.21 .201 1 .21 .91 03. 22.5326. 23. 0880. 24. 0484 .24.8623,25.5747/ 

DATA 1 (ST t 1 .J>, 1*1 . 16). J«=5. 9 I/O. .19.5245. 22. 031 3,23. 4 300.24, 4362. 
125.2221 .26. 1806. 26. 9656, 27.8431 .26.7524,29.51 75.30. 1 789.30.758, 


c 

C GAS MODEL. CONSTANTS 
r 

C»4. 77338* 

Cl ■•00616227 
C2»3. 727440 

C3* «2.*RE0»SCH1/C 1 »+4.*RCM) 

C4> <2.*R£C)/M .+4, »RCH> 

1F»TE'VIP.5T. UO?*,*ALOGIO<P1 + 2600. ) >Gu TO !« 
7 IFIRKO.SO.l .)G0 TO 2 

IF IREO.GT • 1 • )G0 TO 3 • 

C 

C NGN-C I SSOC 1 A T 1 NO GAS MODEL. REQ.LT.l 
C 

X< 1 >■ I2.*C4 t/(C4+2«*C) 
X<2)»I4»*C3)/<C4+2«*C> 
X<3>*l2.-4.«C3-C4>/tC4+2.*C> 
X(3>»<2.*C2»/(C4+2.#C> 

X16)«I4.»C**C1 >/<C4+2.*Cl 
GO TO SO 
C 

C NON-CISSOCtATING GAS model. REO.EO.l 

c 


231 .74.32.55, 23*250* * J6» J 1 * J 7,90* 1 8.90* 1 9*61 . 20 • 1 6 * 20 .62 * 

221 *241 *21 ,tr*9,22.4R7, 22. 912.23.29b. 22. 632, 24, 190 *24 .646*26.031 *0** 
4 16.59*16.6! . 1«.662. 20*401 * 1 6.96*5,21 .632 . 22* 1 56*22. 71 77*23.2772. 
S23.7344.Ji4. ixoe *44.4157 .26.0 1 S *25.4 »4 ,25*6^6 * 0. . 1 1 .280*13.08*14*. 
614.787. 15, 2 W J . 1 9*6. 16.514.17.0724* 17*6301 .16.067.18.4713.18,61 * 

719.3632 * 1 9.81 v r *20. 2043*0*. 20.39. 82. 60 .24.19. 25 .21 .26*007* 26.98. 
627.773*26,66 .29.59. • -‘£.31 .02*31 .61 .32.60.33.45.34.15/ 

CATA ( (ST ( 1 . J ) . ' * 1 , 16 > < J* 10, 12 I/O*. 19.643,21 .2024 ,22.4669, 23.51 56. 
124.300 . 85.234 « 25.98 1 3. 26. 8053. 27. 6588 .28. 3836.29.01 55 .29.5779. 
230.54,21 .36. 22, 06.0.. 21 .62 . 24.31 , 25. 75.26.71 .27.52, 

328. 5805.20.30-'. 30. 269. 31. 250. 32. 036. 32. 71 21 .33.3028.34.3009. 
435.12"0.35.8281 .0. . 1 5.92. 1 7.71 , 16.71 . 19.43. 1 9.979.20.638 .2 l . 1 55 , 
52! .71 3, 22. 2-»2, 22. 727,23. 1126. 23. 4466. 24. 0055. 24. 4654, 24.861 6/ 

C 

C INPUT 
C 

SCO READ (5. 100 1RCH.REQ, TEMP.P 
100 FORMAT 14E16. 8) 

WRITE <6.3001 RCH.REO. TEMP.P 

300 FORMAT! 10X4HRCHsE15.8,4X4HREO>E15.8.4X5HTEMPcE15.8. 1HR. 

1 4X2HP«E 1 5.8, 3HA T M ) 

00 12 J«I,12 
12 X<J)«0. 


2 X < 1 ) s ( £ • * ( l • -2 • *C3 ) )/ ( C+ ( 1 • ~2**C3 ) > 

X (2 )- «2«*C3 >/(C + U • -2.*C3 ) ) 

X(5>*C2/CCM 1 *-2.*C3> ) * 

X(S>=C2**Cc»Cl )/iC+< 1 *-2**C31 ) 

00 TO 00 

3 IF IPCH*E0*0# >00 r C 14 

IF(«r:Q#0 , ‘ •< Cl •+4*^HECH1/(2 **kCH5 ) >GO TO 15 
C 

C MCAI-OISSOCI AT I MG GAS -000L* RtiQ*GT,l 
C 

K0*1C***U .747-<34t5./Ti£«P) 1 
K9~1C*** <?*272- ( 1621C*/TcWP) ) 

AAa< I •-KS>*COC4+2.»C3-I * J**2 

3I3-CC-1 . l-C44.2.#C3i * <2** < l *-2* *C3 ) r.<6?(C4K^.’C3-2. I I 

CC-K8* C2*-C4-4**C3)*2**C3 

X<95- <SG«T<32*»2-4.*AA*CC > )/ <2 • * AA > 

X<2) = C2**C3/€C+C4+2*»C3 *1 . > )-X(9) 

X( 1 )-< C 1 •-2**C1)*X C2 > + C 1 ***C3)»X(9) >/C3 
X ( 4 ) aX < 1 )*X(9)/(<Q*XC2 ) > 

X(5)=C2* C2«*XC2 )+X<V)+XU))/2t 
X<6)=2.*Cl*X(0) 

Pt_! ! *“<9*X(2 )/(X (9)^*2) 

VtfRlTFC6«200)PLIM 


! 


: < 1 =: <694577* /TEMP }-:2- V2v5 ) 

<2=1 0-** C (423o7./TEVP)-6-6325 5 
.<3 = 10. ** ( (47S5C fc/T£y.o ) —6 * 5570 i 
K5=iC«**< (4l6l2*/T£M?;-5i672C} 

<6*1 0*** ( (90047 */T£.V? }-6«.?5C5 ) 

<7=10*** £ <£0437*/TG«P J-6*2v:5} 

IF (RE0*EQ.0* 5 GO 70 4 
I F ( RCH » K'E * G fc ) GO TO 5 
<4=0* .. .. 

GO TC 6 

5 K4 = 10***< <50;73*/7£M?>-7*S765) 

6 A2a=P**3*;<l&;<^(U-£^C3TCi 
A2l«P**2*tKl*C 2 i+C Jt 2* *<2*<4* C 1 *-2**C3i > 
A20=2**F*<2*U .-C3> 

Ai2=P**2*K4*<5* < i .-2**C3-C > 

A1 l=Px . *< l.*C>*2**K4* (1 *-2* vC3* 5 

A 1 0=2 • * ( 1 .-C3> 

AQ3=2 * *P* *2*C*K3*<4 

A02=P*<<4*< 1 .-2**C3+C;*2**X3*<G3+C> J 
A01=l » tC-2.vP*K 4*{ 1 *—2**03 i 
ACC=2 * * ( C3*l * 1 

B22=P**3*< 1 *K4 * C2 * -4 , *C3*-C4 ; 

B2 1 ■P#*2#.(<i * ( 2 • —2» *C3‘-C4 ) + 2 * #<2*<4* < 1 — ? ,*C3) 


! 


. e2C*2**°*K2^ (1 «.-C3) 


212*P**2- <4*<Sv> u i-2v i>C3-C4 ; 


an ap* i <5* < 1 • -C.3-C4 ) +K4* vl-*rS* *C3 > ) 


510=U-C3. 


803*2* *C4*.°*#2#K3*K4 


B02*C4*P* { 2* *rC3+K4 > 


S0l*C4 


C ITERATION FOR CHEMICAL COMPOSITION 


21 Y ( 1 > = % -r I • *SQRT ( C*+C5**£*K3**CvP*K3>-<:»+C5 >✓ (4 *#C*P*K3 1 


I r ( R£Q «. GT 6 * « . wC * C 25 

16 6sCYL}»{Y(ui*ol2+&:i >+31 0 1 -~-*2+4** < Y (L > * < Y >L 5*322+321 )+320l*CY<L5 
1* C Y CL J*<Y(L 5 *303+302 5+301 >*> 

IF<B.LT*0. 5GO TO 20 

Z CL >* CSQR7 CB 5-1** (Y(L>*<YiL 1*512+31 1 i+ciiG > >>< ( Y % - 1* < Y CL ) *322+321 > 
1+6205*2* > 

IF< ( Y CL > *LT« ftOCOC l 5 *At\D« tREQ»LE*l * 5 iGO TO 7 
. F(L,=:Z<L)*{Z{L )*{YCL5*(Y<L> *A22+A2 1 > ~A2C ) - C Y U. 5 * < Y iL } *A 1 2+AI i ) + 


1AI0) >+{Y(L)*<Y(L>*<Y<ui*AG3^A02>+A0: 5+AOO > 


IFC (FvL# *GT *—*0001 ) ® AND* (F*L> aLT» eGCOl 5 5 GO TO S 

CELr?s<2:CL>**2>*C2**A22*YiL)+A2l >+ZCL )* i 2**A 12*Y (L >+A; 1 )+3»*A03* 
1YCL'**2+2**A02*YCL)+A0i- .2**{YCL}*IY 'L >*A22+A21 i +A20 >*Z ) + < Y <L 5 

2*<YCL5*A12+A2 I >*A1C5 )*C CZCL , v*21* (2, *^22*Y*L 5+521 5+Z CL')* <2**©12 
3*YLJ^BU > - C 3**B03*Y CL > **2+2 »*302*Y < L ) +60 1 5 5/12** CYCL >*%YiL)*o22 


4+321 >+B205*ZCL>+<YCL>*<YIL>*8l2+&i: >+310 ) 


YCL+1 1«Y<L)-F<L>/ABSCDELF> 


L«l-+1 
60 TO 16 

20 Ylt_l*.l*YfLl • 

IF1CY1L1* LT.. 00001 ). AND. 1HEQ.LE • 1 • > >60 TO 7 

GO TO 16 

25 Z(1 }*1+1.*SQ«T1 (2.*C4*< U-C31 + 1 »+C >**2+16.*P*K**C4* 11 .-C31*lfc4*l l • 
1-C31+1.+C1 l-12.*C4*U.-C31 + l *+Cl 1/14.*P*K2*1C4*11.-C3>+1 .+C> ) 

26 B» (621*2 1L >**2+81 1*21L)-B0I >** -4**1B20*Z1L>**2+B10*Z1L ) »* 1322 

t*ZlL>**2+B12*ZlL>-802> ‘ 

K«I 

TY CK )* ( — 1 «*SQRT (81 — <621 *Z CL >**2+81 1*ZIL 1— B01 1 >/12»*1822*Zlt 1**2 
1 +B12*Z<I_ 1-602 1 I 

27 0EI.«1B22*TY 1K>**2+B21*TV1K 1+620 >*ZlLl**2+tUl2*TYtK>**2+811*TYtK> 

1 +6 1 0 1 *Z 0. 1 - « B03*TY I K 1 **3+B02*T Y « 1 **2+B0 1 *T Y { K 1 1 

TY1K+1 1«TY1K>— DEL/t— 3»*603*TY |IC >**2+2»* 1822*2 1C >**2+8 12*21 L1-B02 1 
1*TV <K 1+ 1821*Z (L 1**2+B1 1*211. 1-601 1 1 


K=K+1 

GO TO 27 . _ 

26 Y1L1»TY1K> 

•% 

F1L>»Z1L1*1Z11->*1YH->*1Y1I.1*A22+A21 1+A201 + 1 V‘|_)*tYtl.>*Al2+Al 1 > + 
lAIOl 1+1 Y1L 1* l v n -l*lYlL. >*A03+A02 1+A01 1+A001 


|F< <f<L>.LT.. 0001 >.ANO. <F<L>.GT.-»0001 ) >GO TO 6 
rFfZ(L>.LT«.OOOOi >ao TO 7 
IFIL.E0.D60 TO 29 

Z<L.+ I >-Z<L»-FIL»*«Z<L-t t-Z(L))/(F(L-| )“F<L) » 


GO TO 28 


GO To 28 
8 X{7J»YCL> 


IF<REO.EO.O. )G0 TO 22 


XI8 )<Z(U 


Xf I »»P**2*K1*XI7»»XC8)**2 


X«3)aP*X3»X«7l**2 


XC10)»P*K5*X<7>*X<8» 


.XC9l»tC3*«Xll 1+2. *XC3>*XI7>*XC I0>» »✓<<!. -2. *C3»*P*K4#Xt7)*|.-C3> 


X(2I»P*K4*X(7>*X«9) 


X C4 > =P*K2*XI8|**2 . 

X(S>*<C2/2*>*l2»*X«2>+xTl">+2.*Xt3>+X«7>+x'c9>4XCi0>T 


X<6»»2.»CI»XI5| 


GO TO 50 


c dissociating air 
c 

4 A£2*0. 

A£l=0. 


A2C=P*;<6*( 1 *+2**01 ) - 


A * 2 = 0 * 

a: i = p#<7* < : *+Ci > 
A ; c = I •tCI . 


A 0 2 3 0 * 


A02 = P*;<3 


A01=l» .... 

ACG=—I * . . , 

B22=G * 

521=0* 

B20*2**?*i <3 
Bl 2 = 0 « 

511=P*K7*<1 » -C2 ) 
. 510-1 * 

503=0* 

. 302-2 & ■*P*<3*C2 


30 1 =C2 


GO TO 21 
22 X ( 12 >■ 2IL) 


xc i 1 5*?*.<7*X<7)*XC 12 ) 

' X<3)*P*<2*X<7:**2 

X<5l=<C2*<2.*X<3}+X<7>+X<i: 5 5 -X (i 1 ) -X { 1 2 > 5 /2 * 

X<6 >«C1*(2**X<5)tXU 1 >i-X( 12)> 

C , . ’ . , . 

C CALCULATION Or GAS' MIXTURE THERMODYNAMIC :P£k TIES 

C // . _ ’ 

SO DC 9 1-1*16 

' IF(TtW?-TT(I))10 r i^9 
* 9 CONTINUE ’ 

)0 CELH=<TEV.P-TTC ) >/^TT( I 5-TTi 1-1 # , 

DELS* (ALOG < TEMP/T T < I - 1 ) ) . / i ALOG i TT (I >/T7 ;i-l 1 ) ) 

M = 0 • _ 

• • H = 0 • 

”"s*oV’~’ . 

00 'll U-l *12 

’if tx( j;V-P*o« >go tc ;i . .. 

«=M~X<J>*MTCJ> 

K = H *t X £ U } * ( M T \ I — 1‘ * U ) "** £ riT ( l *JJ — HT t : — I ■. J J * * 0 j. '-H ) ^ 

£aS+X<U )* < S7 ( l-l * J )+ CST £ I s J t -ST t 1-1 i Ji i *DELS-ALOG <P*X 
„ 1 1 CONTINUE . . 

H'»(1 *98588/M)»H ' 

' S* Vl *9356S/M>*S 
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The performance parameters are given by the following expressions: 


^ - 0 + t)h - Ss 


e s 


L SP f" 


The variation of specific impulse with Mach number for this idealized 
ramjet is presented in Figure 7« The solid curve was calculated using 
the thermodynamic properties of Reference S. The dashed curve was 
calculated using the computer program for the simplified gas model, 
lhe maximum deviation is 1 per cent. 

Stannary of solid carbon and ammonia formation .- The formation of 
solid carbon extract** useful energy from the combustion gas. If this 
phenomenon occurs during a nozzle expansion process, the useful energy 
absorbed in forming solid carbon is not recoverable and a loss in 
performance results. 

The derivation of the nondissociating (T < T cut _ off ) fuel rich 
gas model required one chemical reaction in order that the formation of 
carbon monoxide (CO) could be Included. The equilibrium expression 
resulting from the required reaction is Independent of pressure. 
Consequently, the composition of the gas in this region is independent 
of pressure. This fact simplifies the treatment of solid 
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carbon formation, since the limiting pressure above which solid carbon 
will form is related to temperature and composition only. 

The computer program for the simplified pis model Indicates the 
formation of solid carbon by an error statement. When the limiting 
pressure is exceeded, the error statement is printed, but the program 
computes the thermodynamic properties as though solid carbon had not 
formed. When the limiting equivalence ratio is exceeded, the error 
statement is printed and no calculations are made. The limiting pressure 
for solid carbon formation is readily extracted from the thermodynamic 
calculations . Figures 8 through 10 present the variation of limiting 
pressure with equivalence ratio and temperature for methane -ai . 

* 0.25), kerosene -air (R^ - 0.5), and benzene-air (R^ ■ l) 
combustion product, respectively. For all cases, as 


R eq ~ >R eqii n 

p lim -*° 

since 

*002 

*eq -1 

p lim 

since 

Xqo -» o 


The formation of ammonia is considered only for the hydrogen-air, 
nondissociating, fuel rich system. The composition of the gas under 
these conditions depends upon equivalence ratio only. Therefore, the 
limiting pressure above which the mole fraction of ammonia '8 greater 
the" some arbitrary amount is related to temperature and composition. 

The computer program for the simplified gas model indicates the 
formation of a monia > 0 .01) by aa error statement. When the 

limiting pressure is exceeded, the error statement is printed, but the 


Limiting pressure, P|j m> o 
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Equivoltnc* ratio, R t q 


Figure 8.- Limiting pressure for solid carbon foimation, methane -air 
combustion products (R^ - 0.25)# 3eq,lia " 







Equivalence rstie,&«q 


Figure 9. 


Malting pressure for solid carton f orma tion, kerosene -air 
co m bu stion products = 0*5)> * 3»0. 





Limiting pressure, P| im> otm 
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Figure 10.- Halting pressure for solid carbon formation, benzene -air 
combustion products (R^ = 1.0), = 2 . 5 . 
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program computes the thermodynamic properties as though ammonia had not 
formed. Figure IX presents the variation of limiting pressure with 
equivalence ratio and temperature for hydrogen-air combustion products. 
For all temperatures, as 


8 

t 

p lim - >#0 

since 

% 2 -0 

«eq-l 

p lim -»• 

since 

% 2 -0 
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Figure 11.- Limiting pressure for ssranla formation (Xjq. • 0.01) 
hydrogen-air combustion products (IL^ * 0). 
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CHAPTER IV 

SUtMARI AND CONCLUSIONS 

A simplified equilibrium hydrocarbon-air combustion gas model for 
convenient use in engine cycle computer programs has been presented. 

The generalized gas model presented Includes the effects of dissociation 
and reduces to the special cases of hydrogen-air combustion products, 
as well as dissociating air. The associated computer program can be 
readily incorporated as a subroutine in a general engine cycle computer 
p rogr am or used alone to generate Mollier diagrams. 

With the exception of pure air, this gas model neglects the 
formation of chemical species containing atomic nitrogen, which allows 
a considerable simplification of the solution for chemical composition. 
The effect of 'tills assumption on the thermodynamic properties of 
stoichiometric kerosene-air and hydrogen-air coLibustion products is 
presented. The importance of this assumption is shown in terms of the 
performance of an Idealized subsonic combustion ramjet. Good agreement 
between the simplified model and more comprehensive treatments is 
obtained in the range of temperatures applicable to hypersonic engine 
cycles. 

The computer program has been used to calculate the limiting 
pressure for solid carbon and ammonia formation in fuel rich gas 
mixtures. The results of this calculation are presented as a function 
of temperature and equivalence ratio for the combustion products of 
various hydrocarbon fuels. 
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